Baryon Stopping as a Probe of Deconfinement Onset in Relativistic Heavy-Ion 

Collisions 
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It is argued that an irregularity in the baryon stopping is a natural consequence of a phase 
transition occurring in the compression stage of a nuclear collision. It is a combined effect of the 
softest point inherent in an equation of state (EoS) with a phase transition and a change in the 
nonequilibrium dynamics from hadronic to partonic transport. Thus, this irregularity is a signal 
from a hot and dense stage of the nuclear collision. In order to illustrate this proposition, calculations 
within the three-fluid model were performed with three different EoS's: a purely hadronic EoS, an 
EoS with a first-order phase transition and that with a smooth crossover transition. It is found that 
predictions within the first-order-transition scenario indeed reveal an a strong irregularity in the 
incident energy dependence of the form of the net-proton rapidity distributions in central collisions. 
This behavior is in contrast to that for the hadronic scenario, where the distribution form gradually 
evolve, displaying no irregularity. The case of the crossover EoS is intermediate. Only a weak 
irregularity takes place. Experimental data also exhibit a trend of similar irregularity, which is 
however based on still preliminary data at energies of 20A GeV and 30A GeV. 

PACS numbers: 25.75.-q, 25.75.Nq, 24.10. Nz 
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I. INTRODUCTION 

Onset of deconfinement in relativistic heavy-oin col- 
lisions is now in focus of theoretical and experimental 
studies of relativistic heavy-ion collisions. This problem 
is one of the main motivations for the currently running 
beam-energy scan [l[ at the Relativistic Heavy-Ion Col- 
lider (RHIC) at Brookhaven National Laboratory (BNL) 
and low-energy-scan program Q at Super Proton Syn- 
chrotron (SPS) of the European Organization for Nuclear 
Research (CERN), as well as newly constructed Facility 
for Antiproton and Ion Research (FAIR) in Darmstadt [3[ 
and the Nuclotron-bascd Ion Collider Facility (NICA) in 
Dubna In this paper I would like to argue that the 
baryon stopping in nuclear collision can be a sensitive 
probe of the deconfinement onset. 

In fact, an irregularity in the incident-energy depen- 
dence of the baryon stopping is very natural if the system 
undergoes a phase transition. Let us start with discus- 
sion of the conventional (i.e. one-fluid) hydrodynamics, 
which is applied to the whole process of the nuclear col- 
lision, e.g. from its compression stage to the expansion 
stage up to freeze-out, like it is done in Refs. [HQ. 

The form of the resulting rapidity distribution of net- 
baryons depends on the spatial form of the produced fire- 
ball. If the fireball is almost spherical, the expansion of 
the fireball is essentially 3-dimcnsional which results in a 
peak at the midrapidity in the rapidity distribution. This 
statement is a theorem that can be proved in few lines. 
If a the fireball is strongly deformed (compressed) in the 
beam direction, i.e. has a form of a disk, its expansion 
is approximately 1-dimensional that produces a dip at 
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the midrapidity, which is confirmed by numerous simu- 
lations, see e.g. Ref. Q- In terms of the fluid mechanics 
this is a consequence of interaction of two rarefaction 
waves propagating from opposite peripheral sides of de- 
caying disk toward its center Q . This speculation is very 
similar to that related to the elliptic flow: a strong elliptic 
flow results from a strongly deformed almond-shaped ini- 
tial fireball with the deformation of the resulting momen- 
tum distribution of particles being inverse to the spatial 
deformation of the initial fireball. The physical mech- 
anism here is precisely the same, only the expansion is 
developed in the transverse direction. 

The next question is how this fireball is formed. This 
is already a matter of dynamics in the early compression 
stage of the nuclear collision. A softest point H charac- 
teristic of EoS's with a phase transition plays an impor- 
tant role in this compression dynamics. At low energies 
the softest point is not reached in the collision process, 
the system remains stiff and therefore the produced fire- 
ball is almost spherical. As a result, the baryon rapidity 
distribution is peaked at the midrapidity. When the inci- 
dent energy gets high enough, the softest-point region of 
the EoS starts to dominate during the compression stage, 
the system weaker resists to the compression and hence 
the resulting fireball becomes more deformed, i.e. more 
of the disk shape. Then its expansion is close to the 1- 
dimensional pattern and, as a result, we have a dip at the 
midrapidity. With energy rise, the stiffness of the EoS (in 
the range relevant to compression stage) grows, the sys- 
tem stats to be more resistant to the compression and 
hence the produced fireball becomes less deformed. The 
expansion of this fireball results in a peak or, at least, to 
a weaker dip at the midrapidity as compared to that at 
the "softest-point" incident energy. With further energy 
rise, the initial kinetic pressure overcomes the stiffness 
of the EoS and makes the produced fireball strongly de- 
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formed again, which in its turn again results in a dip at 
the midrapidity. 

Thus, even without any noncquilibrium, we can expect 
a kind of a "peak-dip-peak-dip" irregularity in the inci- 
dent energy dependence of the form of the net-proton 
rapidity distributions. Noncquilibrium also contributes 
to this irregularity. At a phase transformation 1 the 
hadronic degrees of freedom arc changed to partonic 
ones. In particular, it means a change in cross sections 
which govern the nuclear stopping power. Therefore, this 
change also induces a certain irregularity in the stopping 
power. Only the dip at the midrapidity in ultrarclativis- 
tic nuclear collisions has a different origin in the actual 
case of weak baryon stopping as compared with that in 
the conventional hydrodynamics. It occurs because the 
baryon charges of colliding nuclei traverse through each 
other. 

It is important to emphasize that the "peak-dip-pcak- 
dip" irregularity is a signal from the hot and dense stage 
of the nuclear collision. 

In the present paper this qualitative pattern is illus- 
trated by calculations within a model of the three-fluid 
dynamics (3FD) employing three different equations 
of state (EoS): a purely hadronic EoS [ll| (hadr. EoS), 
which was used in the major part of the 3FD simula- 
tions so far [To|, ll2T - fl5| , and two versions of EoS involv- 
ing the deconfinement phase transition fltl ]. These two 
versions are an EoS with the first-order phase transi- 
tion (2-phasc EoS) and that with a smooth crossover 
transition (crossover EoS). The softest points in these 
EoS's are illustrated in Ref. [13] ■ The hadronic EOS 
possesses no softest point - stiffness of the EoS 
changes monotonously. Preliminary results of simula- 
tions with phase transitions have been already reported 
in Refs. [3 EH. Here I present calculations with thor- 
oughly tuned friction forces in the quark-gluon phase, 
which made it possible to reasonably reproduce a great 
number of observables in the incident energy range 2.7 
GeV < ^snn < 62.4 GcV in terms of the center-of- 
mass energy. Details of this calculations and results on 
a great number of bulk observables (rapidity and trans- 
verse spectra, flow observables and multiplicities) for var- 
ious species and a large number of incident energies, their 
comparison with available data will be reported else- 
where. Here I would like to focus on rapidity distribu- 
tions of net-protons in central collisions of heavy nuclei 
in the AGS-SPS energy range. The distribution of net- 
protons is the best probe of the nuclear stopping in the 
absence of data on net-baryons. 

Concerning the fit of the friction force, few comments 
are in order. The friction in the hadronic phase was es- 
timated in Ref. [2(j. Within hadronic scenario (hadr. 
EoS) we had to enhance this friction in order to repro- 
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duce the baryon stopping at high energies Though 
such a enhancement is admissible in view of uncertainties 
of the estimated friction, the value of the enhancement 
looks too high. Indeed, at ^s^n = 17.3 GeV, i.e. at the 
top SPS energy, this enhancement exceeds the factor of 
2. In scenarios with phase transitions there is no need to 
modify the hadronic friction. This can be considered as 
a theoretical argument in favor of such scenarios. 



II. EQUATIONS OF STATE 

The 2-phase and crossover EoS's still differ even at 
very high densities. The latter means that the crossover 
transition constructed in Ref. fl6| is very smooth. The 
hadronic fraction survives up to very high densities. In 
particular, this is seen from Fig. [TJ the fraction of the 
quark-gluon phase (Wqgp) reaches value of 0.5 only at 
very high energy densities. In this respect, this version 
of the crossover EoS certainly contradicts results of the 
lattice QCD calculations, where a fast crossover, at least 
at zero chemical potential, was found [2l|. Therefore, 
a true EoS is somewhere in between the crossover and 
2-phasc EoS's of Ref. [lj]. 

Figure [T] demonstrates that the onset of the phase tran- 
sition in the calculations happens at top-AGS-low-SPS 
energies. Similarly to Ref. [22(, the figure displays dy- 
namical trajectories of the matter in a central box placed 
around the origin r = (0, 0, 0) in the frame of equal ve- 
locities of colliding nuclei: |a;| < 2 fm, \y\ < 2 fm and 
I z | < 7 cm 2 fm, where j cm is Lorcntz factor associated 
with the initial nuclear motion in the cm. frame. Ini- 
tially, the colliding nuclei arc placed symmetrically with 
respect to the origin r = (0, 0, 0), z is the direction of the 
beam. At a given density n#, the zero-temperature com- 
pressional energy, e(ns,T = 0), provides a lower bound 
on the energy density e, so the accessible region is cor- 
respondingly limited. In the case of the crossover EoS 
only the region of the mixed phase between Wqgp = 0.1 
and Wqgp — 0.5 is displayed, since in fact the mixed 
phase occupies the whole (e-ns) region. The s-ns repre- 
sentation is chosen because these densities are dynamical 
quantities and, therefore, are suitable to compare calcu- 
lations with different EoS's. 

Only expansion stages of the evolution are displayed, 
where the matter in the box is already thermally equi- 
librated. Evolution proceeds from the top point of the 
trajectory downwards. Subtraction of the tojvtib term 
is taken for the sake of suitable representation of the 
plot . The size of the box was chosen to be large enough 
that the amount of matter in it can be representative to 
conclude on the onset of the phase transition and to be 
small enough to consider the matter in it homoge- 
neous medium. Nevertheless, the matter in the box still 
amounts to a minor part of the total matter of colliding 
nuclei. Therefore, only the minor part of the total matter 
undergoes the phase transition at 10A GeV energy. 

As seen, the phase transition starts at the top AGS en- 
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FIG. 1: Dynamical trajectories of the matter in the central box of the colliding nuclei (4fmx4fmX7 cm 4fm), where y C m is the 
Lorentz factor associated with the initial nuclear motion in the cm. frame, for central (b = 2 fm) collisions of Au+Au at 4A 
and 10 A GeV energies and Pb+Pb at 20 A GeV. The trajectories are plotted in terms of baryon density (jib) and the energy 
density minus ns multiplied by the nucleon mass (e — mjvns). Only expansion stages of the evolution are displayed. The 
trajectories are presented for two EoS's: 2-phase EoS (left panel) and crossover EoS (right panel). Symbols on the trajectories 
indicate the time rate of the evolution: time span between marks is 1 fm/c. For the 2-phase EoS (left panel) the shadowed 
"mixed phase" region is located between the borders, where the QGP phase start to raise (Wqgp — 0) and becomes completely 
formed (Wqgp = 1). For the crossover EoS (right panel) the corresponding borders correspond to values of the QGP fraction 
Wqgp = 0.1 and Wqgp = 0.5. Inaccessible region is restricted by e(ns, T = 0) — mjvns from above. 



crgies in both cases. It gets practically completed at low 
SPS energies in the case of the case of the 2-phase EoS. 
In the crossover scenario it lasts till very high incident 
energies. The trajectories for two different EoS's are nev- 
ertheless very similar at displayed energies. Apparently, 
it happens because the friction in the quark-gluon phase 
was selected in such a way that both scenarios reasonably 
reproduce available data at high energies. 



III. PROTON AND NET-PROTON RAPIDITY 
DISTRIBUTIONS 

A direct measure of the baryon stopping is the net- 
baryon (i.e. baryons-minus-antibaryons) rapidity distri- 
bution. However, since experimental information on neu- 
trons is unavailable, we have to rely on net-proton (i.e. 
proton-minus-antiproton) data. Presently there exist ex- 
perimental data on proton (or n et-p roton) rapidity spec- 
tra at AGS [HHp and SPS HMof energies. These data 
were analyzed within various models [!, H, [hJ [H, HH-HH . 

Figure [5] presents calculated rapidity distributions of 
protons (for AGS energies) and net-protons (for SPS en- 
ergies) and their comparison with available data. Notice 
that difference between protons and net-protons, as well 
as a contribution of weak decays to these yilds are negligi- 
ble at the AGS energies, see compilation of experimental 
data in Ref . [4(| . Contribution of weak decays of strange 



hypcrons into proton yield was disregarded in accordance 
with measurement conditions of the NA49 collaboration. 
Correspondence between the fraction of the total cross 
section related to a data set and a mean value of the im- 
pact parameter was read off from the paper [HJ in case of 
NA49 data. For Au+Au collisions it was approximately 
estimated proceeding from geometrical considerations. 

As seen from Fig. [2j at lower AGS energies all EoS's 
predict the same results, since at these energies only 
hadronic parts of all EoS's are relevant, see Fig. [T] Re- 
sults of the 2-phase EoS start to differ from those of the 
hadr. and crossover EoS's beginning from 6A GeV: the 
2-ph.-EoS distributions reveal a dip at the midrapidity. 
This dip contradicts the available experimental data and 
is very robust: variation of the friction in a wide range 
does not remove this dip. Therefore, it is a direct con- 
sequence of the onset of the first-order phase transition, 
which starts precisely at these energies in the 2-ph.-EoS 
scenario, see Fig. [1] This dip survives even in one-fluid 
calculations involving the lst-order phase transition 
in spite of immediate baryon stopping inherent in the 
one-fluid model. In 3FD calculations, this dip changes 
into midrapidity peak at higher energies (30A GeV and 
40 A GeV). With further energy rise (E lab > 40 A GeV) 
the midrapidity peak again transforms into a dip, see also 
Fig. [21 The latter dip is already a normal behavior which 
takes place at arbitrary high energies. 

The experimental distributions exhibit a qualitatively 
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FIG. 2: Rapidity spectra of protons (for AGS energies, see left panel) and net-protons (p — p) (for SPS energies, see right panel) 
from central collisions of Au+Au (AGS) and Pb+Pb (SPS). Experimental data are from collaborations E895 [HI, E877 [2^ |. 
E917 [25|], E866 [2(|, and NA49 [27l43ll |. The percentage shows the fraction of the total reaction cross section, corresponding 
to experimental selection of central events. Feedback of weak decays into p and p yields is disregarded. 



similar behavior as that in the 2-phase-EoS scenario: af- 
ter a plateau or a shallow midrapidity dip at the energy 
of 8A and 10A GeV, a well pronounced peak at 20A GeV 
is again observed. However, quantitatively the 2-phasc- 
EoS results certainly disagree with data at 8A GeV, 10A 
GeV and 40A GeV energies. They also disagree with data 
20A GeV and 30A GeV, which however have preliminary 
status, and hence it is too early to draw any conclusions 
from comparison with them. The 2-phase-EoS behavior 
is in contrast with that for the hadronic-EoS scenario, 
where the form of distribution in central collisions grad- 
ually evolve from peak at the midrapidity to a dip. The 
case of the crossover EoS is intermediate: only a shallow 
dip occurs at 10A and 15A GeV while at 20A GeV the 
distribution looks like a plateau. 

Predictions of different scenarios diverge to the largest 
extent in the energy region 8A GeV < E[ a b < 40A GeV. 
Unfortunately data at 20A and 30A GeV still have a 
preliminary status and disagree with any considered sce- 
nario. Updated experimental results at energies 20A and 
30A GeV are badly needed to pin down the preferable 
EoS and to check the hint to the zigzag behavior of the 



type "peak-dip-peak-dip" in the net-proton rapidity dis- 
tributions. 



IV. ANALYSIS OF "PEAK-DIP-PEAK-DIP" 
IRREGULARITY 

In order to quantify the above-discussed "peak-dip- 
peak-dip" irregularity, it is useful to make use of the 
method proposed in Ref. [l8[ . For this purpose the data 
on net-proton rapidity distributions are fitted by a simple 
formula 

= a (exp{-(l/w s )cosh(y - y cm - y s )} 

+ cxp{-(l/w s )cosh(y-y cm + y s )}) (1) 

where a, y s and w s arc parameters of the fit. The form 
(fT]) is a sum of two thermal sources shifted by ±y s from 
the midrapidity. The width w s of the sources can be 
interpreted as w s = (temperature (/(transverse mass), if 
we assume that collective velocities in the sources have 
no spread with respect to the source rapidities ±y s . The 
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Vs NN [GeV] Vs NN [GeV] Vs NN [GeV] 

FIG. 3: Midrapidity reduced curvature [see. Eq. (J2j)] of the (net)proton rapidity spectrum as a function of the center-of-mass 
energy of colliding nuclei as deduced from experimental data and predicted by 3FD calculations with different EoS's: the 
hadronic EoS (hadr. EoS) [ll|] (left panel), the EoS involving a first-order phase transition (2-ph. EoS, middle panel) and the 
EoS with a crossover transition (crossover EoS, right panel) into the quark-gluon phase Upper bounds of the shaded areas 
correspond to fits confined in the region of \y — y cm \/ycm < 0.7, lower bounds — \y — y C m\/ycm < 0.5. 
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FIG. 4: The same as in Fig. [3]but for midrapidity value of the 
(net)proton rapidity spectrum scaled by y cm - The thin long- 
dashed line corresponds to the hadr. -EoS calculation without 
fragment production, i.e. without coalescence. 



parameters of the two sources are identical (up to the 
sign of y s ) because only collisions of identical nuclei are 
considered. 

The above fit has been done by the least-squares 
method and applied to both available data and results 
of calculations. The fit was performed in the rapidity 
range \y — y cm \ly C m < 0.7. The choice of this range is 
dictated by the data. As a rule, the data are available 
in this rapidity range, sometimes the data range is even 



more narrow (80 A GeV and new data at 15SA GeV |3l1|). 
I put the above restriction in order to treat different data 
in approximately the same rapidity range. Another rea- 
son for this cut is that the rapidity range should not be 
too wide in order to exclude contribution of cold specta- 
tors. The fit in the rapidity range \y — y C m\/ycm < 0.5 
has been also done in order to estimate uncertainty of 
the fit parameters associated with the choice of fit range. 

A useful quantity, which characterizes the shape of the 
rapidity distribution, is a reduced curvature of the spec- 
trum at the midrapidity defined as follows 



dN\ 



= (y cm /w s ) 2 (sinh 2 y s - w s coshy s ) 



(2) 



The factor 1/ {y C mdN/dy) m is introduced in order to 
get rid of overall normalization of the spectrum. The 
second part of Eq. @ presents this curvature in terms 
of parameters of fit |T]). The reduced curvature, C y , and 
the midrapidity value, (y C mdN/dy) y _ y ^ m , are two inde- 
pendent quantities quantifying the the spectrum in the 
midrapidity range. Excitation functions of these quan- 
tities deduced both from experimental data and from 
results of the 3FD calculations with different EoS's are 
displayed in Figs. [3] and 2] Notice that a maximum in 
y cm {dN/dy) cm at ^snn = 4.7 GeV happens only be- 
cause the light fragment production becomes negligible 
above this energy. The 3FD calculation without coales- 
cence (i.e. without the light fragment production) reveals 
a monotonous decrease of y C m(dN/dy) cm beginning from 
the lowest energy considered here. 
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FIG. 5: The same as in Fig. [4] but for midrapidity values of the net-proton (left panel) and proton (right pan el) rapidity 
spectrum in conventional re pre sentation, i.e. without scaling. Experimental data are from collaborations E895 |23(|. E877 [24| . 
E917 [U, E866 Sl> NA49 [HH and STAR g2. 



To evaluate errors of C y values deduced from data, I es- 
timated the errors produced by the least-squares method, 
as well as performed fits in different the rapidity ranges: 
\y - VcmMVcm < 0.5 and \y - y cm \/y cm < 0.7, where it 
is appropriate. Problems were met in fitting the data at 
80 A GeV [H and the new data at 158 A GeV [Hj]. These 
data do not go beyond the side maxima in the rapidity 
distributions. This results in large uncertainty in the pa- 
rameters. In particular, because of this problem I keep 
the old data at 158^4 GeV [27J in the analysis. The er- 
ror bars present largest uncertainties among mentioned 
above. The upper errors of 158A-GeV 80A-GeV points 
results from the uncertainty of the narrow rapidity range. 
The uncertainty associated with the choice of the rapid- 
ity range turned out to be a dominant one in the case 
computed data. Therefore, in Fig. [3] results for the com- 
puted spectra are presented by shaded areas with borders 
corresponding to the fit ranges \y — y C m\/ycm < 0.7 and 
\y — ycm\/ycm < 0.5. In Fig. 0]the midrapidity values of 
the rapidity spectra were taken directly from experimen- 
tal data and calculated results. Therefore, only experi- 
mental error bars are displayed there. 

Since experimental data at AGS and RHIC energies 
were taken from Au+Au collisions while at SPS the 
Pb+Pb collisions were studied, the calculations were per- 
formed respectively for Au+Au (6=2 fm) and Pb+Pb 
(b = 2.4 fm) central collisions. In fact, at the same in- 
cident energy the computed results for Pb+Pb collisions 
at b = 2.4 fm are very close to those for Au+Au at b = 2 
fm. Therefore, the corresponding irregularity of the en- 
ergy dependence of the fit parameters is negligible. 

The irregularity in data is distinctly seen here as a 
zigzag irregularity in the energy dependence of C y . Of 



course, this is only a hint to irregularity since this zigzag 
is formed only due to preliminary data of the NA49 col- 
laboration. A remarkable observation is that the C y exci- 
tation function in the first-order-transition scenario man- 
ifests qualitatively the same zigzag irregularity (left panel 
of Fig. [3]) as that in the data fit, while the hadronic 
scenario produces purely monotonous behavior. The 
crossover EoS represents a very smooth phase transition, 
as mentioned above. Therefore, it is not surprising that 
it produces only a weak wiggle in C y . 

This zigzag irregularity of the first-order-transition sce- 
nario is also reflected in the midrapidity values of the 
(net)proton rapidity spectrum (Fig. 0]). In the conven- 
tional representation of the data without multiplying by 
y cm , the irregularity of the (dN/dy) cm data is hardly 
visible (Fig. [5]). Moreover, the difference between pre- 
dictions resulted from different EoS's is also scarcely dis- 
cernible. Thus, the scaled representation of Fig. [4] serves 
as a zoom revealing differences. In Fig. [5] not only net- 
protons but also proton midrapidity values are displayed 
in a wider energy range. However, results for top energy 
y/sNN = 62.4 GeV are very approximate, since more ac- 
curate computation requires unreasonably high memory 
and CPU time. As seen, a visible difference between 
net-protons and proton data, as well between predictions 
of hadronic EoS and EoS's with phase transitions starts 
only at RHIC energies. Similar (but based on differ- 
ent, double-gaussian fit) analysis was performed in Rcf. 
[4lj . It found no irregularities in excitation functions of 
the corresponding parameters, in particular, because the 
analysis of Ref. [4JJ was restricted to energies Ei a b > 
20A GeV. 

As it was pointed out in the introduction, the "peak- 
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dip-peak-dip" irregularity is very natural in a system un- 
dergoing a phase transition. First, it is associated with 
the softest point of a EoS. Therefore, it is not surprising 
that the irregularity is weaker in the crossover scenario 
than in the first-order-transition one. Indeed, the softest 
points in the crossover EoS is less pronounced than in the 
first-order-transition one (l7j . There is no softest point 
in the hadronic EoS and hence there is no irregularity. 

The second reason of this irregularity is a change in 
the noncquilibrium regime. The 3FD model takes into 
account the leading noncquilibrium of the nuclear colli- 
sion associated with a finite stopping power of the nuclear 
matter. It simulates the finite stopping power by means 
of friction between three fluids. Naturally, this friction 
changes when a phase transition happens. In the case 
of the crossover scenario this change in the friction is 
very smooth. Therefore, it does not contribute to the 
irregularity. At the same time this change in the friction 
enhances the irregularity in the first-order-transition sce- 
nario. As it was demonstrated in Ref. [l8|, if the same 
friction is used in both phases, the reduced curvature cal- 
culated with the 2-phasc EoS reveals a wiggle behavior in 
C v but with considerably smaller amplitude as compared 
with zigzag in actual calculations with different frictions 
in different phases. These different frictions appear quite 
naturally in the 3FD model. The hadronic friction was 
estimated in Ref. [2(| and works well at lower AGS en- 
ergies. Therefore, there are no reasons to modify it. The 
partonic friction, while not microscopically estimated, is 
fitted to reproduce data at high incident energies. This 
is a reason to believe that it is a proper choice. 



V. CONCLUSION 

An irregularity in the baryon stopping is a natural con- 
sequence of a phase transition occurring in the compres- 
sion stage of a nuclear collision. It is a combined effect of 
the softest point of a EoS and a change in the noncquilib- 
rium regime from hadronic to partonic one. It is impor- 
tant to emphasize that this irregularity is a signal from 
the hot and dense stage of the nuclear collision. 

In order to illustrate this fact, calculations within the 
3FD model were performed with three different equa- 
tions of state, i.e. with a purely hadronic EoS [ll|, and 
also with two versions of EoS involving the deconfine- 
ment phase transition jTij : an EoS with the first-order 
phase transition and that with a smooth crossover transi- 
tion. The crossover transition constructed in Ref. [l6| is 
very smooth. The hadronic fraction survives up to very 
high energy densities. In this respect, this version of the 
crossover EoS certainly contradicts results of the lattice 
QCD calculations, where a fast crossover, at least at zero 
chemical potential, was found [2l[ . Therefore, a true EoS 
is somewhere in between the crossover and 2-phasc EoS's 
of Ref. [H. 

It is found that predictions within the first-order- 
transition scenario reveal a "peak-dip-peak-dip" irregu- 



larity in the incident energy dependence of the form of 
the net-proton rapidity distributions in central collisions. 
At low energies, rapidity distributions have a peak at the 
midrapidity. With the incident energy rise it transforms 
into a dip, then again into a peak, and with further en- 
ergy rise the midrapidity peak again change into a dip, 
which already survives up to arbitrary high energies. The 
behavior the type of "peak-dip-pcak-dip" in central colli- 
sions within the 2-phasc- EoS scenario is very robust with 
respect to variation of the model parameters in a wide 
range. This behavior is in contrast with that for the 
hadronic-EoS scenario, where the distribution form grad- 
ually evolve from peak at the midrapidity to a dip. The 
case of the crossover EoS is intermediate. Only a weak 
wiggle of the type of "pcak-dip-peak-dip" takes place. 

Experimental data also revel a trend of the "peak- 
dip-peak-dip" irregularity in the energy range 8A GeV 
< Eiab < 40A GeV, which qualitatively similar to that 
in the first-order-transition scenario while quantitatively 
differ. However, the this experimental trend is based on 
preliminary data at energies of 20A GeV and 30^4 GeV. 
In general, predictions of different scenarios diverge to 
the largest extent in the energy region 8A GeV < Ei a b < 
40A GeV. Therefore, updated experimental results at en- 
ergies 20A and 30A GeV are badly needed to pin down 
the preferable EoS and to check the hint to the zigzag 
behavior of the type "pcak-dip-peak-dip" in net-proton 
rapidity distributions. Moreover, it would be highly de- 
sirable if new data in this energy range are taken within 
the same experimental setup and at the same central- 
ity selection. Hopefully such data will come from new 
accelerators FAIR at GSI and NICA at Dubna. 

In order to quantify the observed "peak-dip-peak-dip" 
irregularity, the a naly sis of the distribution shape pro- 
posed in Refs. (l8l.ll9| was applied. This method is based 
on calculation of the reduced curvature of the spectrum 
at the midrapidity C y . In terns of C y the irregularity in 
data is distinctly seen as a zigzag irregularity in the en- 
ergy dependence of C y . The energy location of this zigzag 
anomaly coincides with the previously observed anoma- 
lies for other hadron-production properties at the low 
SPS energies 0, |45[. A remarkable observation is that 
the C y energy dependence in the first-order-transition 
scenario manifests qualitatively the same zigzag irregu- 
larity, as that in the data fit, though quantitatively does 
not reproduce the data fit. The hadronic scenario pro- 
duces purely monotonous behaviour. The crossover EoS 
represents a very smooth phase transition, as mentioned 
above. Therefore, it is not surprising that it results in 
only a weak wiggle in C y . 

Here the distribution of net-baryons was discussed be- 
cause net-baryons have advantage of being confined by 
the baryon number conservation. In general, an analy- 
sis of the form of particle spectra looks very promising. 
For instance, in Ref. [43| it was proposed to analyze 
the incident-energy dependence of the width of pion ra- 
pidity distributions. This width can be associated with 
the sound velocity in the dense stage of the reactions. 



8 



It is found that the sound velocity has a local minimum 
(indicating a softest point in the EoS) in the range of 
tJsnn = 4-9 GeV, which coincides with the range of the 
"peak-dip-peak-dip" irregularity discussed in this paper 
and that of the previously observed anomalies for other 
hadron-production properties 0, EBJ . 
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